New potential super-incompressible phase of ReN2 
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Abstract 

The structural, elastic, and electronic properties of ReN2 are investigated by first-principles 
calculations with density functional theory. The obtained orthorhombic Pbcn structure is energet- 
ically the most stable structure at ambient pressure. ReN2 is a metallic, superincompressible solid 
and presents a rather elastic anisotropy. The estimated Debye temperature and hardness are 735 
K and 17.1 GPa, respectively. Its estimated hardness is comparative to that of SisN4. 

PACS numbers: 61.50.Ks, 62.20.-x, 71.15.Mb, 71.20.-b, 71.20.Be 
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I. INTRODUCTION 



? or these reasons, most recently, three transition metal dinitrides, PtN2, 1] 
J] have been successfully synthesized under extreme conditions of pressure 



Transition metal nitrides are of great technological and fundamental importance because 
of their strength and durability as well as their useful optical, electronic, and magnetic 
properties. [l| By incorporating B, C, N, or O atoms into the interstitial sites of transition 
metals with highly valence electron density (VED), the compounds of transition metal with 
highly hardness and strong incompressibility can be synthesized by means of high pressure 

IrN 2 ,fl and OsN 2 , 

and temperature. Very large bulk moduli observed in experiments make them be candidates 
of superhard materials, which could be used in cutting tools and wear-resistant coatings. 
Meanwhile, many theoretical investigations have been performed to explore their structures, 
for that the crystal structure is an important prerequisite of understanding their physical 
properties. The search for the structural forms of PtN 2 , OsN 2 and IrN 2 with the lowest 
energy and mechanical stability usually considers hypothetic structures of orthorhombic, 
tetragonal, hexagonal and cubic crystal lattice. The structures of PtN 2 and OsN 2 are indexed 
as cubic (pyrite structure, Pa3) [l|, |5j and orthorhombic (marcasite structure, Pnnm) . 
While it is very difficult to determine the structure of IrN 2 in experiment because of the large 
intensity ratio (1:100) of X-ray diffraction peaks between nitride and pure metal. The first- 
principles calculations have shown that IrN 2 is monoclinic (CoSb 2 structure, P2 1 /c).[8( As a 
neighbor of noble transition metal Os, the diboride of Re has been successfully synthesized,^, 
9[ however the dinitride of Re to date has not been reported experimentally. Theoretically, 
Zhao et al suggested that ReN 2 with P4 2 /mnm structure is thermodynamically stable at 



ambient conditions and up to 76 GPa. lOj Here, we readdress structural and elastic properties 
of ReN 2 and find that a new phase with Pbcn symmetry is more stable than that with 
P^/mnm symmetry. 

II. COMPUTATIONAL DETAILS 



All calculations are performed by the CASTEP codepjj using ab initio pseudopotentials 
based on DFT with the exchange-correlation functional of Ceperley and Alder as parameter- 
ized by Perdw and Zunger (LDA-CAPZ).[l3] All the possible structures concerned are opti- 
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mized by the BFGS algorithm (proposed by Broyden, Fletcher, Goldfarb, and Shannon). [13l| 
which provides a fast way of finding the lowest energy structure and supports cell optimiza- 
tion in the CASTEP code. The interaction between the ions and the valence electrons is 
described by using Vanderbilt's supersoft pseudopotential[14]. We used a plane wave cutoff 
of 310 eV and a Brillouin zone sampling grid spacing of 2tcx 0.04 A -1 . Pseudoatomic calcu- 
lations are performed for N 2s 2 2p 3 and Re 5s 2 5p 6 5d 5 6s 2 . In the geometrical optimization, all 
forces on atoms are converged to less than 0.002 eV/A, all the stress components are less than 
0.02 GPa, and the tolerance in self-consistent field (SCF) calculation is 5.0 xlO -7 eV/atom. 
Relaxation of the internal degrees of freedom is performed at each unit cell compression or 
expansion. 

The standard method is used to calculate the elastic constants, in which the second 
derivatives of the internal energy of a crystal are determined as a function of properly chosen 
lattice distortions describing strains. In practice, by applying small elastic strains, the elastic 
constants are obtained from the change of energy or stress. Elastic moduli are given by using 
the elastic constants according to the Voigt-Reuss-Hill (VRH) approximation, [jjj] The Hill 
average, in general, is selected as the estimation of bulk modulus and shear modulus. 



III. RESULTS AND DISCUSSIONS 



A. Structural property 

To search for the most stable structure of ReN 2 , we consider monoclinic, orthorhombic, 
tetrahedral, hexagonal, and cubic structures with a Re:N stoichiometry of 1:2, including 
these initial structures of AuTe 2 (C2/m, No. 12), CoSb 2 (P2i/c, No. 14), GeS 2 (Fdd2, No. 
43), CaCl 2 (Pnnm, No. 58), OsB 2 (Pmmn, No. 59), /?-Re0 2 (Pbcn, No. 60), Pbca (8c 
Wycknoff sites, No. 61), Pnma (4c Wycknoff sites, No. 62), ReSi 2 (Immm, No. 71), SiS 2 
(Ibam, No. 72), Cu 2 Sb (PA/nmm, No. 129), rutile (P4 2 /mnm, No. 136), MoSi 2 (IA/mmm, 
No. 139), P3ml (la and 2d Wycknoff sites, No 164), CrSi 2 (P6 2 22, No. 180), Fe 2 P (P62m, 
No. 189), A1B 2 {PQ/mmm, No. 191), ReB 2 (P6 3 /mmc, No. 194), pyrite (Pa3, No. 205), 
fluorite (Fm3m, No. 225), and NiTi 2 (Fd3m, No. 227). Combining total energy with 
elastic constants analysis, four competing structures, i.e., Pbcn, P4 2 /mnm, Pmmn, and 



Fm3m, satisfy the mechanical stability critetion, 
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171 ] in which the latter three phases 
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have been discussed by Zhao et al.|10j The equilibrium volume Vq and bulk modulus Bo 
are determined by fitting the total energy as a function of volume to the 3rd-order Birch 
Murnaghan equation of state (EOS). [3] Our calculated equilibrium parameters (see Table I) 
show that the rutile (P4 2 /mnm) and the Pbcn phases have larger equilibrium volumes than 
the fluorite (Fm3m) phase. From the calculated relative enthalpy AH, one can see that the 
Pbcn phase is the most stable structure at zero pressure. In the optimized Pbcn structure, 
four rhenium atoms occupy 4c Wycknoff sites (y=0.1191) and eight N atoms hold 8d sites 
(x=0.2663, y=0.3526, z=0. 0871). The nearest distance between Re (N) and N atoms is 
1.955 (2.644) A. The orthorhombic Pbcn structure possesses a structure characterized by 
zigzag chains of Re atoms propagating along the c axis of the unit cell with bond length and 
bond angle of 2.794 A and 122.4°, respectively. In order to determine the structural phase 
transition between competing structures, the enthalpy H of each phase of ReN2 is computed. 
The enthalpy difference per formula unit (f.u.) as a function of pressure is plotted in Fig. 
1. We do not observe structural phase transition from Pbcn phase to others up to 100 GPa. 
It is interesting to notice that Pbcn structure has been observed in experiment being the 
ground state of ReOo. 

VED is defined as the total number of valence electron divided by volume per formula 
unit, which is an important factor when searching the superhard materials. The valence 
electron shell of Re is 6s 2 and 5d 5 , and that of N is 2s 2 p 3 so that the total number of valence 
electrons is 17 for ReN2- VED obtained in four phases (Table I) is slightly lower than that 
of Re metal (0.4761 electrons/A 3 ) due to the adding of N atoms. High VED is good to 
resist the fracture, which contributes to high bulk pressure of ReN2 (Table I), predicting 
that ReN 2 could be considered as a candidate of a superhard material. 



B. Elastic property 

The elastic constants of Pbcn structure obtained from our calculations (cn=365 GPa, 
c 22 =553 GPa, c 33 =610 GPa, c 44 =138 GPa, c 55 =230 GPa, c 66 =83 GPa, ci 2 =307 GPa, 
Ci3=232 GPa, C23=242 GPa) satisfy the mechanical stability requirement well, implying 
that the orthorhombic Pbcn structure is mechanically stable. As is concluded above, ReN 2 
has a very high bulk modulus, and therefore can be regarded as a candidate of a superhard 
material. However, a high bulk modulus by itself does not directly lead to high hardness. In 
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general, both high bulk modulus and high shear modulus are indispensable to high hardness. 
Therefore, the shear modulus must be taken into account when searching for new superhard 
materials. The bulk modulus represents the resistance to fracture, while the shear modulus 
represents the resistance to plastic deformation. Moreover, two other factors are important 
for technological and engineering applications: Young's modulus E and Poisson's ratio v. 
Young's modulus, defined as the ratio between stress and strain, is used to provide a measure 
of stiffness of the solid. The larger the value of E, the stiffer the material. The Young's 
modulus and Posson's ratio for an isotropic material are given by E = ^^ G , v = ^b+g) > 
respectively, [151] where B and G represent Hill average bulk modulus and shear modulus. 
The calculated bulk modulus B, shear modulus G, and Young's modulus E of ReN 2 with 
Pbcn structure are 334 GPa, 127 GPa, and 337 GPa, respectively, of which the bulk modulus 
of ReN2 is slightly smaller than that of Re metal (360 GPa) . It is worthy to note that the 
bulk modulus obtained from the elastic constants agrees well with the one obtained through 
the fit to the 3rd-order Birch Murnaghan EOS (Bo), providing a consistent estimation of the 
compressibility for ReN2. Possion's ratio v reflects the stability of a crystal against shear. 
This ratio can formally take values between -1 and 0.5, which corresponds respectively to 
the lower limit where the material does not change its shape, or to the upper limit when 
the volume remains unchanged. Poisson's ratio of Pbcn structure is 0.3319, which is larger 
than 0.3, indicating that there exists small volume change during elastic deformation. Also 
Possion's ratio provides more information on the characteristics of the bonding forces. It has 
been proved that z/=0.25 is the lower limit for central-force solids and 0.5 is the upper limit, 
which corresponds to infinite elastic anisotropy. The high value v for ReN2 means that the 
central interatomic forces are central. [2^] In order to predict the brittle and ductile behavior 
of solids, Pugh introduced the ratio of the bulk modulus to shear modulus of polycrystalline 
phases. 2l|] A high (low) B/G value is associated with ductility (brittleness). The critical 
value which separates ductile and brittle materials is about 1.75. The significant difference 
obtained between bulk and shear modulus indicates the higher ductility of ReN2- 

Now, we turn to discuss the elastic anisotropy of ReN2, which is estimated by the universal 
anisotropy index A = 5^- + ^-— 6. [22] The discrepancy of A from zero determines the extent 
of single crystal anisotropy and accounts for both the shear and the bulk contributions. The 
obtained A value is 1.252, indicating that ReN2 has a large elastic anisotropy. Additionally, 
fractional axis compression as a function of pressure are plotted in Fig. 2. The curves 
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presented in this figure clearly exhibit that: a bigger difference in fractional axis compression 
values exists for Pbcn phase, implicating that there is a evident structure deformation with 
increasing pressure; there is a transition point for compression along b- and c-axis with 
increasing pressure, that is, fe-axis is more difficult to be compressed than c-axis at 38 
GPa above. This abnormality for compression along different axis can foretell that there 
is a structural phase transition at about 40 GPa. Further analysis of the change of the 
structural parameters (atomic coordination, bond length, et al.) shows that from 40 to 45 
GPa fractional coordinations x and z in 8d (x, y, z) sites (N atoms) change much, yielding 
the great change of the nearest distance between N atoms (cIn-n)- In detail, at 40 GPa, x 
and z are 0.2184 and 0.0752, respectively. However, x and z reduce to 0.1385 and 0.0349 at 
45 GPa, respectively. Correspondingly, the d N _ N reduces from 2.359 A at 40 GPa to 1.432 
A at 45 GPa. 

Further, the Debye temperature Op which correlates with many physical properties of 
materials, such as specific heat, elastic constants, and melting temperature is given. The 

1 /3 

Debye temperature Qd can be calculated by the equation 0£> = f [f^ (^m 1 )] [f (<§) 3//2 



|( g+4G/3 ) 3//2 ]~ 1//3 ' 23( | where h is Plank's constant, k is Boltzmann's constant, Na is Avo- 
gadro's number, p is density, M is the molecular weight and n is the number of atoms in the 
molecule. The calculated Gd value of 735 K in Pbcn is comparative to the values of 866.4 K 



in ReB2[24j, 850 K in Re02 19j and 691 K in OsN2. 6| Employing the correlation between 
the shear modulus and Vickers hardness reported by Teter [25] for a wide variety of hard 
materials, the indentation hardness of ReN 2 can be estimated to be approximately 17.1GPa 
, which is comparable to that of SiaN4 (21±3 GPa). 2ij] It is well known that the hardness of 
superhard materials should be higher than 40 GPa, indicating that ReN 2 is not a superhard 
material. Jhi et al proposed that the hardness is determined by elastic constant C44 
rather than shear modulus. In our case, we note that there is little difference between the 
C44 value of elastic constant and shear modulus. Accordingly, both Teter estimation and 
Jhi's viewpoint of hardness actually achieve a consistent result for ReN 2 . 

C. Electronic property 

In order to understand the nature of elastic property, the electronic density of states is 
also analyzed. The total and partial density of states (DOS) with respect to Fermi level 
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is shown in Fig. 3. Finite DOS at Fermi level N(Ep) of 5.035 states/eV shows that ReN2 
is metallic. The electrons from Re-5ci and the N-2p states both contribute much to the 
DOS near the Fermi level. That is, the electronic structure of ReN 2 is governed by strong 
hybridization between the Re-5<i and N-2p states, imaging that there is a covalent bond 
in between Re and N atoms. The existence of covalent bonding contributes to large bulk 
modulus, i.e., strong incompressibility. However, we have not found any pseudogap that 
observed in ReB 2 ,Q, Q OsB 2 



and OsN 2 ,[6f implying that covalent bonding in between 
Re and N atoms is not so strong. Besides, the calculated net charges on N and Re atom 
are -0.57 e and 1.15 e, which indicates that the chemical bonding between Re and N have 
some characteristics of ionicity. The occurrence of ionic bonding results in the smaller shear 
modulus, which explains the reason of the not so large hardness in ReN 2 . 



IV. CONCLUSION 



In conclusion, the most stable structure of ReN 2 is explored based on first-principles 
methods in the framework of density functional theory within local density approximation. 
It results that orthorhombic Pbcn structure is the most stable structure at ambient pressure. 
ReN 2 is a metallic, super-incompressible solid and presents large elastic anisotropy. The 
estimated hardness is comparative to that of Si3N 4 . The structural determination and the 
discussion of the elastic and electronic properties may be of use in the synthesis and high 
pressure study of ReN 2 . 



V. ACKNOWLEDGEMENT 



This work was supported by the special Funds for Major State Basic Research Project 
of China(973) under grant no. 2007CB925004, 863 Project, Knowledge Innovation Program 
of Chinese Academy of Sciences, and Director Grants of CASHIPS. Part of the calcula- 
tions were performed in Center for Computational Science of CASHIPS and the Shanghai 



7 



Supercomputer Center. 



[1] J. C. Crowhurst, A. F. Goncharov, B. Sadigh, C. L. Evans, P. G. Morrall, J. L. Ferreira, and 

A. G. Nelson, Science 311 (2006) 1275. 
[2] R. B. Kaner, J. J. Gilman and S. H. Tolbert. Science 308 (2006) 1268. 

[3] H-y. Chung, M. B. Weinberger, J. B. Levine, A. Kavner, J-M. Yang, S. H. Tolbert, and R. B. 

Kaner, Science 316 (2007) 436. 
[4] A. F. Young, C. Sanloup, E. Gregoryanz, S. Scandolo, R. J. Hemley, and H. K. Mao, Phys. 

Rev. Lett. 96 (2006) 155501. 
[5] A. F. Young, J. A. Montoya, C. Sanloup, M. Lazzeri, E. Gregoryanz, and S. Scandolo, Phys. 

Rev. B 73 (2006) 153102. 
[6] Z. Wu, X. Hao, X. Liu, and J. Meng, Phys. Rev. B 75 (2007) 054115. 

[7] Z. W. Chen, X. J. Guo, Z. Y. Liu, M. Z. Ma, Q. Jing, G. Li, X. Y. Zhang, L. X. Li, Q. Wang, 

Y. J. Tian, and R. P. Liu, Phys. Rev. B 75 (2007) 054103. 
[8] R. Yu, Q. Zhan, and L-C. De Jonghe, Angew. Chem. Int. Ed. 46 (2007) 1136. 
[9] S. La Placa and B. Post, Acta Crystallogr. 15 (1962) 97. 
[10] E. Zhao and Z. Wu, Comput. Mat. Sci. 44 (2008) 531. 

[11] M. D. Segall, P. J. D. Lindan, M. J. Probert, C. J. Pickard, P. J. Hasnip, S. J. Clark, M. C. 

Payne, J. Phys.: Cond. Matt. 14 (2002) 2717. 
[12] D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45 566 (1980); J. P. Perdew and Y. Wang, 

Phys. Rev. B 45 (1992) 13244 . 
[13] B. G. Pfrommer, M. Cote, S. G. Louie, and M. L. Cohen, J. Comput. Phys. 131 (1997) 133. 
[14] D. Vanderbilt, Phys. Rev. B 41 (1990) R7892. 
[15] R. Hill, Proc. Phys. Soc. London 65 (1952) 349. 

[16] O. Beckstein, J. E. Klepeis, G. L. W. Hart, and O. Pankratov, Phys. Rev. B 63 (2001) 134112. 
[17] M. Born and K. Huang, Dynamical Theory of Crystal Lattices, Clarebdom, Oxford, 1956. 
[18] J. -P. Poirier, Introduction to the physics of the Earth's Interior, Cambridge University Press, 

Cambridge, 2000. 
[19] Y. L. Li and Z. Zeng, Chin. Phys. Lett 25 (2008) 4086. 

[20] P. Ravindran, L. Fast, P. A. Korzhavyi, B. Johansson, J. Wills, and O. Eriksson, J. Appl. 

8 



Phys. 84 (1998) 4891. 
[21] S. F. Pugh, Philos. Mag. 45 (1954) 823. 

[22] S. I. Ranganathan and M. Ostoja-Starzewski, Phys. Rev. Lett. 101 (2008) 055504 . 
[23] O. L. Anderson, J. Phys. Chem. Solids 24 (1963) 909. 

[24] X. Hao, Y. Xu, Z. Wu, D. Zhou, X. Liu, X. Cao, and J. Meng, Phys. Rev. B 74 (2006) 224112. 
[25] D. M. Teter, MRS. Bull. 23 (1998) 22. 

[26] S-H. Jhi, J. Ihm, S. G. Louie, and M. L. Cohen, Nature 399 (1999) 132. 
[27] Y. L. Li, G. H. Zhong, and Z. Zeng, accepted by Chin. Phys. B (2009) 

[28] H. Y. Gou, L. Hou, J. W. Zhang, H. Li, G. F. Sun, and F. M. Gao, Appl. Phys. Lett. 88 
(2006) 221904. 



9 



TABLE CAPTIONS 



TABLE I: Equilibrium lattice parameters, V (A 3 ), a (A), 6(A), c (A), density p (g/cm 3 ), 
valence-electron density p e (electrons/A 3 ), bulk modulus B (GPa), and enthalpy difference 
per chemical formula unit (f.u.) AH (eV). V is of per chemical f.u. 
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FIGURE CAPTIONS 

Fig.l: Enthalpy difference per formula unit as a function of pressure. P4 2 /mnm is taken 
as a reference point. 

Fig. 2: Fractional axis compression as a function of pressure for Pbcn structure. 

Fig. 3: Total and partial density of states for Pbcn structure. Vertical dotted line indicates 
the Fermi level. 
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TABLE I: 



Li et al. 



Pbcn P '42 jmnm Fm3m Pmmn Reference 



V 


31.251 


31.693 


27.712 


26.779 


This work 






31.7 


27.7 


26.7 


[10] 


a 


4.518 


4.769 


4.801 


8.088 


This work 






4.775 


4.803 


8.083 


[10] 


b 


5.657 






2.380 


This work 










2.385 


[10] 


c 


4.902 


2.789 




2.782 


This work 






2.78 




2.77 


[10] 


P 


11.385 


11.217 


12.838 


12.406 


This work 


Pe 


0.416 


0.410 


0.469 


0.453 


This work 


B 


331.3 


337.8 


381.9 


365.2 


This work 


AH 





0.015 


1.457 


1.663 


This work 
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FIG. 1: Li et al. 
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FIG. 2: Li et al. 
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FIG. 3: Li et al. 
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